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Abstract 



In this paper we introduce the space-fractional Poisson process whose state probabilities p^{t), t > 0, 
a e (0, 1], are governed by the equations (d/dt)Ps:(t) = — A°(l — B)p^(t), where (1 —By is the fractional 
difference operator found in the study of time series analysis. We explicitly obtain the distributions p"it), 
the probability generating functions G^iu, t), which are also expressed as distributions of the minimum of 
^ I i.i.d. uniform random variables. The comparison with the time-fractional Poisson process is investigated 

and finally, we arrive at the more general space-time fractional Poisson process of which we give the 
explicit distribution. 
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1 Introduction 

Fractional Poisson processes studied so far have been obtained either b y considering renewal processes with 



Fractional Poisson processes studied so rar nave been oDtained eitner d v considering renewal 
inte rtimes between events represented by Mittag-Leffler distributions I^ Mainardi et alj,l2004l. 



Beghin and Orsinghei . 



I2OO9I or by replacing the time derivative in the equations governing the state probabilities with the frac- 
tional derivative in the sense of Caputo. 
In this paper we introduce a space-fractional Poisson process by means of the fractional difference operator 

A« = (l- Br, aG(0,l], (1.1) 



which often appears in the study of long memory time series [iTsavl. boosll . 

The operator dl.iP implies a dependence of the state probabilities p^it) from all probabilities p'^^t), j < k. 

For a = 1 we recover the classical homogeneous Poisson process and the state probabilities P]^lt) depend 

onlyonpfc_i(t). 



For the space-fractional Poisson process we obtain the following distribution: 

p«(t) = Pr{iV«(t) = fc} = ^2^ -±—^—-J- k>0. (1.2) 

fc! -^ r! rfar + 1 — fc) 

The distribution of the space-fractional Poisson process can be compared with that of the time -fractional 
Poisson process N^.{t), t > 0, v e (0, 1] : 

Pr{iV,(t) = fc} = ^— ^^ ^ ^ fc>0. (1.3) 

fc! ^ r! r(v(fc-|-r)-|- 1) 

For a = V = 1, from (I1.2D and (I1.3D . we immediately arrive at the classical distribution of the homogeneous 

Poisson process. 

The space-fractional Poisson process has probability generating function 

G^iu, f) = Eu'^"^'^ = e-^"(i-")"^ \u\ < 1, (1.4) 

and can be compared with its time-fractional counterpart 

,Giu,t) = Eu''^^'^=E,.i-Xa-u)tn, \u\<l, (1.5) 

where 

is the one-parameter Mittag-Leffler function. 

We show below that the probability generating function of the space-time fractional Poisson process reads 

,G„(u,t) = £,,(-A"(l-urt^'), |u|<l, (1.7) 

and its distribution has the form 

(-1)'' ~ (-A^t")'- r(ar + l) 

Pk^'^ = -T^T.7^^^7F^^T^^' ^ > 0, ae (0,1], v e(0,l]. (1.8) 

" fc! ^ r(vr -I- Ij r[ar -I- 1 — fcj 

We also show that the space-fractional Poisson process JV"(t) can be regarded as a homogeneous Poisson 
process iV(t), subordinated to a positively skewed stable process S'^^t) with Laplace transform 

lEg -ijs«(t) ^ g-t.«^ z > 0, t > 0. (1.9) 

In other words, we have the following equality in distribution 

iV"(f) = JV(S°(f)). (1.10) 

The representation (II.IOD is similar to the following representation of the time -fractional Poisson process: 

JV,(t) = JV (72,(0), (1.11) 

where T2,,(t), t > 0, is a process whos e one-dimensional distribution is obtained by folding the solution to 

the time-fractional diffusion equation fBeghin and Orsingher, 2009] . 

Finally we can note that the probability generating function ( 11.71 1, for all u e (0, 1), can be represented as 

,.G„(u,0 = Pr| min xy">l-ul, (1.12) 

where the Xj^s are i.i.d. uniformly distributed random variables. 



2 Construction of the space-fractional Poisson process 

In this section we describe the construction of an alternative fractional generalisation of the classical ho- 
mogeneous Poisson process. First, let us recall some basic properties. Let us call 



P,(t) = Pr{Ar(t) = fc} = e 



fc! ' 



t > 0, A > 0, fc > 0, 



(2.1) 



the state probabilities of the classical homogeneous Poisson process Nf^t), t > 0, of parameter A > 0. It is 
well-known that the probabilities Pi^^t), k>0, solve the Cauchy problems 



^''^ ^ ll, fc = 0. 



(2.2) 



Starting from ( 12. 2D . some time -frac tional general is ations of the homogeneous Poisson process have bee n 
introduced in the literature (see e.g. lLaskinl [|2003ll . lMainardi et all [|2004l . lBeghin and Orsinghen [|2009fl ). 



These works are based on the substitution of the integer-order derivative operator appearing in (12.21) with 
a fract ional-order derivative operator, such as t he Riemann-Liouville f raction al derivative (as in jLaskin 
II2OO3I] ) or the Caputo fractional derivative (as in lBeghin and Orsingherl [|2009l ]). In this paper instead, we 
generalise the integer-order space-difference operator as follows. First, we rewrite equation (12. 2D as 



^p,(t) = -A(l-B)p,(t), 



(2.3) 



where B is the so-called backward shift operator and is such that B{p^{t)) = pj:_i(t) and B''(p^(t)) = 
B''~^{B{pj^{t))) = P(._r(t). The fractional difference operator A" = (1 — B)" has been widely used in 
time series analysis for constructing processes displaying long memory, such as the autoregressive fraction- 
ally integrated moving average process (ARFIM A). F'or mor e information on long memory processes and 
fractional differentiation, the reader can consult Tsa^ [I2OO5II . page 89. 
Formula (12. 3D can now be easily generalised by writing 



^^p^M = -A°(l - B)X(0, a e (0, 1], 



(2.4) 



where p^^t), k > 0, t > 0, represents the state probabilities of a space-fractional homogeneous Poisson 
process N'^it), t >0, i.e. 



p°(t) = Pr{Af"(t) = fc}, k>0. 
In turn, we have that ( I2.4D can also be written as 

■ ip^.W = -A« Y!:-=o TJ^}^,i--Lypt-rM, a e (0, 1], 



(2.5) 



P,"(0) ■ 



0, fc > 0, 

1, fc = 0. 



(2.6) 



Note that, in (12.61 ) we considered that p'^^t) = 0, ; e Z . Equation (12.61) can also be written as 
^^p«(t) = -AX(t) 



ip^M = -AX(t) + aAX-i(t) - ^Pk-2M + ■■■+ (-i)'c+i «(«-i)-4^-'c+i) pg(t)^ 

(2.7) 



By applying the reflection property of the gamma function r(z)r(l — z) = 7i/sin(7rz) for z = r — a, we 
have also that 

f^p^it) = -X'^plit) + aA>«_,(t) + i^^ i:'=2 B(a + 1, r - a)p«_^(t), 

(2.8) 



,'^-^{:: ::: 



where the sum is considered equal to zero when fc < 2 and B(x,y) is the beta function. From ( 12.71 ) and 
(12. 8D we see that for a = 1 we retrieve equation (12. 6D of the homogeneous Poisson process. 
An example of process whose state probabilities PkM, depend on all Pjit), j < k, is the iterated Poisson 
process N^t) = JViCJVjCt)), where N^ (t) and Nr,( t ) are i ndependent homogeneous Poisson processes. The 



process JV(f) is analysed in lOrsingher and Politol 11201 111 . 



Theorem 2.1. Let N'^{t), t > 0, be a space-fractional homogeneous Poisson process of parameter A > and 
let G^(^u, f) = Eu'' '~'^\ \u\ < 1, a e (0, 1], be its probability generating function. The Cauchy problem satisfied 
by G^(u, t) is 

|^G>,t) = -A«G„(u,t)(l-u)«, |iz|<l, 
\g„(u,0)=1, 

with solution 

G„(u,t) = e-^"^(i-")", |u|<l, (2.10) 

that is, the probability generating function of a discrete stable distribution (see ' Devrovei U99^1 . page 349). 
Proof. Starting from ( 12. 6D . we have that 

d CO 00 f — 1)'' 

— Gjiz, t) = -A"r(a + l)XiXi'^' irr^i T^"-^^^ ^^.11) 

ot „, r!r(a + l — r) "^ "^ 

= _A«r(a + 1) 2 2 ,^, \' M ) 
^—^j—ir\Y{a-\-l — r) 



°° f— 1) 

= -rna+l)G,[u,t)Y, 



^^^r\na + l-r) 
= -A«G„(u,t)(l-u)«, 

thus obtaining formula ( I2.9D . It immediately follows that 

G>,0 = e-^"^^i-"^", |ii|<l. (2.12) 

D 



Remark 2.1. Note that, for a = 1^ formula (12.101 ) reduces to the probability generating function of the 
classical homogeneous Poisson process. Furthermore, from (I2.10D we have that E[JV"(t)]^ = oo, j e N, 
a e (0, 1). 



Remark 2.2. LetX^, k = l,..., be i.i.d. Uniform[Q, 1] random variables, then 



Gaiu, t) = e 



-?.'Ki-ur_p^l min Xl'''>l-u 

0<k<N[t) "^ 



ue(0,l), (2.13) 



where Af(tX t > 0, is a classical homogeneous Poisson process of rate X"^ with the assumption that min(Xj^) = 
1 when Nit) = 0. 

Theorem 2.2. The discrete stable state probabilities of a space-fractional homogeneous Poisson process N^^t), 
t > 0, can be written as 



p«(t)=Pr{iV«(t) = fc} 



(2.14) 



(-1)*^ ^ i-X^ty r(ar + l) 
-^ r! rfar + l- 



fc! 
fc! 



li)! 



a, a) 
(l-/c,a) 



r(ar + l-fc) 



t > 0, fc > 0, 



where ^i/'jC^) '■s the generalised Wright function (see \Kilbas et al.\ \200di] . page 56, formula (1.11.14)). 
Proof. By expanding the probability generating function (I2.10I) we have that 

-A°t(l-ii)° 



G^iu, t) = e" 



(2.15) 



r=0 

^ (-A"t)'- ^ (-u)"'r(ar + 1) 
^ r! Alm!r(ar + 1- 



r=0 ■ m=0 



;^o" ^'- h ''■ r(ar + l- 



r(ar + 1 — m) 

ar - 
r(ar + 1 — m) 



From this, formula ( I2.14D immediately follows. 



D 



Remark 2.3. Note that the discrete stable distribution ( 12.141 1 (which for a = 1 reduces to the Poisson distri- 
bution) can also be written as 



^'^ '-I ^ r! rfar + l- 



fc! 
fc! 



r=0 

JC /-OO 00 



r=0 



r(c 
!r(c 



fc) 



(2.16) 



A^r!r(ar + l- 



fc) 



Theorem 2.3. Let S^(^t), t > 0, y e (0, 1), be a j-stable subordinator, that is a positively skewed stable process 
such that 



Ee 



-zsr^t) ^ g-tzi- 



z > 0, t > 0, 



(2.17) 



and with transition function q^is, t). For a space-fractional Poisson process N^^t), t > 0, a e (0, 1]^ with rate 
A > 0, the following representation holds in distribution: 

JV«(S''(f)) = iV"''(t). (2.18) 

Proof. In order to prove the representation (I2.18D it is sufficient to observe that 

G^^u,s)qJs,t)ds= I e-^''^(i-"5"q,(s,t)ds = e-^°'''fi-")°'. (2.19) 



r 





Remark 2.4. JVote that, when a = l, formula (I2.18D reduces to 



m\ ^—f: r\ rfar + 1 — ml 



Hence 



and therefore 



'^ ^ m! ^ (m-1)! 

m=k m=k 

-^ ml -^-J 



, m! -^^ m! 

fc m=Jc-l 



(fc-1)! 



D 



NiSKt))=NKt), (2.20) 

and this reveals a second possible way of constructing a space-fractional Poisson process. 

Consider now the first-passage time at k of the space-fractional Poisson process 

T^(t) = inf{t:JV"(t) = fc}, fc>0. (2.21) 

Since Pr{T^ < t] = Pr{Af°(t) > fc}, we have that 

~ (-l)-" ~ (-A«t)'" r(ar-M) 
Pr{T«<f} = y^ ^y^ r^ -^- (2.22) 



Pr{T° e ds} = y ^^ — ^y^^ ^ —■ (2.23) 

^ '' ^ m! ^ (r-1)! r(ar -I- 1 - m) 

m=k r=l V y ^ J 

Note that for a = 1 we obtain the classical Erlang process. First we have 

Pr{Tl<t} = Y,-^J^T '- (2.24) 

■^ m! ^—' (r — mj! 

m=k r=m ^ ■' 

^ m! ' 



(2.25) 



Remark 2.5. Note that, with a little effort, fractionality can be introduced also in time. In this case, for 
example, the fractional differential equation governing the probability generating function is 



j^.G^iu, t) = -A«,,G„(u, t)(l - uT, \u\ < 1, V G (0, 1], a e (0, 1], 



(2.26) 



where d^' /dt^' is the Caputo fractional derivative operator (see [Kilbas et al.\ \200dl ). By means of Laplace 
transforms, some simple manipulations lead to 

„G,iu,t) = E,i-reil-uri \u\<l, (2.27) 

where E^,{x) is the Mittag-Leffler function ^Kilbas et all \200dil . In turn, by expanding the above probability 
generating function we have that 



a,v , (-1)'' ^ (-A«t')'- r(ar + l) 



fc! ^r[vr + l)r[ar + l-k)' 



fc>0, ae(0,l], ve(0,l]. (2.28) 



When a = 1 these probabilities easily reduce to those of a fractional Poisson process (see \Beghin and Orsingher 






fc! ^r(vr + l)r(r-fc + l) 
fc! A'TYvr 



(2.29) 



r=k 



T{vr + l){r-k)\ 



S-ir 



r-Jr\ i^n- 



kj r(vr + l) 



Moreover, let Xj^, k = 1,..., be Ltd. Uniform[0,l] random variables and A/,,(t), t > 0, be a homogeneous 



time-fractional Poisson process of rate A" with the assumption that min(X, ) = 1 when iV^.(t) = 0. The 



probability generating function ,,G„(u, t), for u e (0, 1), can be written as 
,G„(u,t) = £,(-A°r(l-u)°) 

r=0 



(2.30) 



r(vr + l) 



~ (A«t^)'- ' .fr\ , 

Z^-1)'"-Fr— irS^-l^ ^ [l-(l-li)«]* 
^ r(vr + l)^ ^kj 

Y.i^-n-urfJ^i-iy-'' ^- — - 



k=0 



r=k 



kJ r(vr + l) 
2 [Pr [xl'" >l-u) ] 'Pr {JV„(t) = fc} 



k=0 

-- Pr \ min x}'" > 1 - u !■ . 

[0<J:<N,,(t) '^ 



Formula (I2.30t shows that the contribution of the space-fractionality affects only the uniform random variables 
X, while the time-fractionality only the driving process N^XO- 
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